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The Higgs boson production can be affected by several anomalous couplings, e.g. Ct and Cg 
anomalous couplings. Precise measurement oi gg ^ h production yields two degenerate parameter 
spaces of Ct and Cg; one parameter space exhibits the SM limit while the other does not. Such a 
degeneracy could be resolved by Higgs boson pair production. In this work we adapt the strategy 
suggested by the ATLAS collaboration to explore the potential of distinguishing the degeneracy at 
the 14 TeV LHC. If the Ct anomalous coupling is induced only by the operator HQlHIh, then 
the non-SM-like band could be excluded with an integrated luminosity of ~ 235 fb~^. Making use 
of the fact that the Higgs boson pair is mainly produced through an s-wave scattering, we propose 
an analytical function to describe the fraction of signal events surviving a series of experimental 
cuts for a given invariant mass of Higgs boson pair. The function is model independent and can be 
applied to estimate the discovery potential of various NP models. 


Introduction: In the establishment of the Standard 
Model (SM), thousands of experiments have been per¬ 
formed to measure model parameters and check the con¬ 
sistency of the theory. In the same spirit, it is critical 
to measure all the properties of the recently discovered 
Higgs boson as precisely as possible to test the SM and 
probe New Physics (NP) beyond the SM. The Higgs cou¬ 
plings have been constrained at the LHC Run-1 [ma¬ 
in the SM the Higgs boson is predominantly produced 
through the gluon fusion process which can be affected 
by either hti or hgg anomalous couplings. The anoma¬ 
lous couplings are described by the effective Lagrangian 
C = -^cthit + . Figure [^displays the 

parameter space of Ct and Cg allowed by the measure¬ 
ment. Two degenerate parameter spaces arise from the 
interference of c* and Cg contributions. For example, in 
the heavy top-quark limit, 

ctnp( 55 ^ /i) ^ crsuigg h)x \cg + ctp. (1) 

The upper band in Fig. [^corresponds to Ct + Cg ^ -fl, 
which has a SM limit oi Ct —t 1 and Cg —t 0. We name 
it as “SM-like” band. The lower band corresponds to 
Ci + Cg ^ —1, which does not exhibit the SM limit. In 
the lower band the NP contributions should cancel the 
SM contribution out and leads to a residual contribution 
of minus SM value. Though often ignored, it is possible 
in principle We name the lower band as the “faked-no- 
new-physics” (FNNP) band since further improvement in 


^ For example, additional colored SU(2)l singlet scalars (Si) 
could generate a large negative Cg in the FNNP band. The scalar 
Si’s interact with the SM Higgs boson via —kiSfSiHlH with 
H the SM Higgs boson doublet. Integrating out heavy Si’s in¬ 
side the gg h triangle loop yields Cg ~ Y■■ TfS'd kjv’^ /(4m|, ) 
where T{Si) is the Dynkin index of Si for the corresponding 
representation under SU{3)c- 



Ct 

FIG. 1. The bounds at 95% confidence level (C.L.) in the 
plane of Ct and Cg from the ATLAS [1] and CMS [2] collabo¬ 
rations. The red box denotes the SM (ct = l,Cg = 0). 


the measurement of single Higgs-boson productions can¬ 
not distinguish it from the upper SM-like band. The 
Higgs boson pair production gg —>■ hh, which is highly 
correlated with the gg ^ h process, can be used to dis¬ 
criminate the SM-like and FNNP bands [SH^. In this 
work, we focus our attentions on the ATLAS constraints 
and explore the potential of the 14 TeV LHC to exclude 
the FNNP band in the Higgs boson pair production. 
Higgs boson pair production: The Higgs boson pair 
production is usually considered as the best channel to 
measure the Higgs trilinear coupling in the SM [5H13]- It 
is also sensitive to various NP models pUlbSj . In this 
work we adapt the effective Lagrangian approach to de¬ 
scribe the unknown NP effects. After the electroweak 
symmetry breaking the effective Lagrangian related to 
the double Higgs production is nissHn] 


Ch = - —CthtRh - ^c^hhtBh^ 

Zv V 






127n; 


247n;2 


Cr + h.c. , (2) 


where v = 246 GeV is the vaccum expectation value, 
with gs the strong coupling strength, mt is 






2 


the top-quark mass and mh is the Higgs boson mass. In 
the SM, C 3 = Ct = 1 and C 2 h = Cg = 0. The squared am¬ 
plitude of gg —)■ hh averaging over the gluon polarizations 
and colors is [53] 


IMP = 


2567r2u4 


iml 


rC3 


Ct^A +0^9 


+ CtF\2 + -Cg 


+ IctCnl 


2c2/tFA 


(3) 


where Fa = Fa(s, t, ), Fn = Fn(s, f, m^) and 

Ga = Gn(s, t, ruj) are the form factors | 21 | with s 
and t the canonical Mandelstam variables. Gn corre¬ 
sponds to the d-wave component which is negligible |46j . 

In order to compare (j{gg —>■ hh) with the SM predic¬ 
tion, we define a ratio Rhh as 


Rhh 


<^{99 -» hh) 

cr^^(gg —>■ hh) 


(4) 


Figure [^displays the contours of Rhh = 1,5 and 10 in the 
plane of anomalous couplings for four benchmark choices: 


(a) C 3 = 1 , C 2 h = 0 ; (b) C 3 = 1 , C 2 h = 3(ct - l)/ 2 ; 

(c) C 3 = ct = 1; (d) Ct = 1, C 2 h = 0. (5) 


In the case (b) we assume that both ct and C 2 h anomalous 
couplings are induced by the operator HQ^HIr and 
thus exhibit the relation. The production cross section 
of the Higgs boson pair is enhanced in the FNNP band 
in all the four choices: 


(a) 1.44 < Rhh < 92.89; ( 6 ) 11.12 < Rhh < 46.58; 
(c) 3.46 < Rhh < 88 . 86 ; (d) 9.11 < Rhh < 12.64. 


Collider simulation: Next we perform a detailed Monte 
Carlo simulation to estimate the needed integrated lumi¬ 
nosity for probing or excluding the FNNP band at the 
14 TeV LHC. As a concrete example, we examine the 
hh —>■ 6677 channel, which has been studied by the AT¬ 
LAS collaboration m- MadGraph5 m is used to gen¬ 
erate the signal events at the parton-level with CT14 m 
and MSTW2008 [75] parton distribution function (PDF). 
Following the ATLAS study the signal events must 
contain two 6 -tagged jets and two isolated photons which 
satisfy the kinematic cuts as follows: 

^teading b ^ ^ 25 GeV, \r]^\ < 2.5, 

pj. > 30 GeV, \r]^\ < 1.37 or 1.52 < < 2.37, 

AFq < ^Rhh^-yy < 2.0 , ARhy > AFq, AFq = 0.4 , 

lOOGeV < uihi < 150GeV, > llOGeV, 

123GeV < < 128GeV, > llOGeV. 

In order to mimic the imperfect detector effects, we smear 
the final state parton momenta by a Gaussian distri¬ 
bution as suggested in Ref. [75]. The b identification 
strongly depends on p^ and 77 ^. We fit the 6 -tagging ef¬ 
ficiency given in the ATLAS Technique Report m and 
obtain the following function of the 6 -tagging efficiency: 

e, (PT, V) = 0.135 tanh tanh 

X [3 + g-(l’7l-VPT’/iooo)"/i-6jg-hl='PT/iooo^ 

where pr is in the unit of GeV. Figure [^displays our 6 - 
tagging efficiency as a function of p^ and 77 ^, which agrees 
well with the ATLAS study. 






FIG. 2. The contours of Rhh = 1, 5 and 10 in the plane of 
anomalous couplings at the 14 TeV LHC. The figure indices 
correspond to the benchmark choices of anomalous couplings 
in Eq. [^ 
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FIG. 3. The Flagging efficiency as a function of and 

The photon energy resolution and identification are 
crucial to trigger and reconstruct the signal events. In 
the simulation we adapt the photon energy resolution 
and identification efficiency given by the ATLAS collab¬ 
oration [77]. The photon energy resolution is given as 
follows: 

CT (GeV) = 0.3 © 0.10 x y/E{GeV) © 0.010 x F(GeV), 
for |? 7 | < 1.37, 

CT (GeV) = 0.3 © 0.15 x A/F(GeV) © 0.015 x F(GeV), 
for 1.52 < \r]\ < 2.37. (7) 
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The identification efficiency, which is sensitive to p^, is 
given by 


(pt) = 0.76 - 1.98 exp (- , (8) 


Note that the identification rate is less than 80% even for 
an energetic photon. 

The key of collider simulation is to know the so-called 
cut efficiency, i.e. the fraction of signal events passing 
the kinematic cuts. To understand the cut efficiencies 
of different values of anomalous couplings, one has to 
repeat the collider simulation which include all the kine¬ 
matic cuts, imperfect detector resolutions and particle 
identifications, etc. However, it is very time consuming 
in practice. Inspired by the scalar feature of Higgs boson, 
we propose an analytical function to describe the fraction 
of signal events passing through the kinematic cuts. The 
function depends upon the invariant mass of the Higgs 
boson pair (mhh) and is not sensitive to those anomalous 
couplings or specific NP model. The advantage of our 
method is that the cut efficiency of the hh signal events 
can be easily obtained from the convolution of the dif¬ 
ferential cross section of da/drrihh and the cut efficiency 
function. The method is explained below. 

The scattering of gg —)■ hh is dominated by the s- 
wave contribution. Owing to the scalar feature of the 
Higgs boson, there is no spin correlations among the ini¬ 
tial state and final state particles. Therefore, the px and 
r] distributions of the 6-jets and photons depend mainly 
upon mhh- The differential cross section of the gg —> hh 
process before any cut is 


da 

dnihh 


mhh 


n {mhh.P^r) 


^fg/p {xi,gf) 




mhh 


n {mhh, fj^r) s {mhh, s, Pf) . 


(9) 


where ry and i) is the rapidity of one of the Higgs bosons 
in the laboratory frame and center of mass (c.m.) frame, 
respectively. y/S is the collision energy of the hadron 
collider, 'H{mhh, Pr) is the hard scattering cross section 
with pr the renormalization scale, fg/p is the gluon PDF 
with pf the factorization scale. As argued above, the 


cut efficiency depends on the configuration of the Higgs 
bosons which is described by mhh and ry, and the differ¬ 
ential cross section after cuts is 


dac 


dmhh 


[ drhhh'^^^'H {fhhh, Pr) f — 

' fg/p ^ fg/p {^l,hf) J 


dfj 

dr] 


e{mhh,mhh,xi,r]), ( 10 ) 

where fhhh is introduced to describe the finite energy 
smearing effect and e stands for the cut acceptance. De¬ 
fine 


E {rhhh,S,pf) = / dm, 


hh 


dxi 

js 


'xfg/p{xi,Pf) I dp . 

Imhh,?7,a;i 

xe{mhh,mhh,xi,p), (11) 


fg/p 

dfj 

dp 


'-hh 


\XiS 


hf 


then 

^CUt “ 


drhhh^^^hi{rhhh,hr)^{mhh,S,pf). ( 12 ) 


We introduce a differential cut efficiency function 

( 13 ) 


which depends on V~S, parton distribution functions and 
kinematic cuts. Again, we emphasize that the A-function 
is model independent. If the gg —)■ hh scattering in NP 
models is also dominated by the s-wave scattering, its 
production rate after cuts is 

7 da . 

dmhhp -<8 A 

dmhh 

The A-function could be derived analytically for a given 
mhh- It is hard to account for the detector effects analyt¬ 
ically, however. In this work we first obtain the cut effi¬ 
ciency of signal events from the Monte Carlo simulation 
with all the detector effects and then fit the efficiency 
with the following functions as suggested by analytical 
calculations: 


{mhh,S,pf). (14) 



A{mhh) = < 


1 


1 - 


m\h (1 ~ cos Ai?o) — 8 {mh — Smi) 


\(1 


— cos AR, 


i) {r 


1 - 




T,h 


fhh -4,{mh- 6m2) 


^hh ~ 4 {mh - dmi) 
mhh \ 


VsJ 


1 + /3c 


mhh \ I ‘2mhh '^ 


Vs 


Vs ;j’ 


^ (*) 
r^hh > m\l^, 


329.3GeV < mhh < m^^\, 


mhh < 329.3GeV. 


(15) 
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FIG. 4. (a): Cut efficiency as a function of rrihh at the 14TeV LHC; (b): comparison between our method (solid curve) and 
the ATLAS simulation (dashed line, see Fig. 8 in [Z2]); (c) comparison of 68 % C.L. exclusion contours in the plane of Cg and 
ct obtained in our method (solid curve) and the one in [1] (dashed curve) for the choice (b), i.e. C 2 h = 3(ct — l)/2. 


The 7 c, /3a,h,c and dmi ^2 parameters reflect the imperfect 
detector effects, ci _2 are the normalization parameters 
and is the turning point of two fitting functions. 

Figure Qa) displays the cut effldency as a function of 
rrihh from the Monte Carlo simulation result (circle). We 
note that the efficiency depends mainly on the cut on 
and AR. The large cuts require the Higgs 

boson’s energy must be more than + p^ ~ 167 GeV 
such that rrihh > 330 GeV. Therefore, the cut efficiency 
is zero for mhh < 330 GeV. The efficiency increases with 
rrihh because the pr's of Higgs boson decay products also 
increase with rrihh such that more signal events pass the 
Pt cuts. On the other hand, decreases with rnth 

and reaches Ai?o = 0.4 at rrihh = i^hh- i^hh > "^^hh 
the cut efficiency decreases with rrihh because Ai?^^ is 
likely to be smaller than 0.4 such that most of the signal 
events fail the AR cut. We fit the Monte Carlo data with 
the ^-function in Eg. |15| and obtain those fitting param¬ 
eters which are shown in Table |Tj We note that the fitting 
parameters are not sensitive to the PDF sets. Both CT14 
and MSTW2008 PDF sets yield similar parameters. For 
comparison we also plot the fitting function in Fig.|^a); 
see the solid curve. 


TABLE 1. Fitting parameters (the second and fifth rows) of 
A{mhh) and the corresponding uncertainties (the third and 
sixth rows) for two PDF sets: (top) CT14 [74] and (bottom) 
MSTW2008 [75]. 


rrihh ,0 

Cl mi 7c 

C2 5m2 pa Pb Pc 

1.257TeV 

1.1378 50GeV 1.675 

11.02 2.5GeV 1.13 1.48 4.88 

- 

0.0037 - 0.002 

1.20 - 0.02 0.02 0.03 

rrihh ,0 

Cl mi 7c 

C2 Sm2 Pa Pb Pc 

1.254TeV 

1.1419 50GeV 1.677 

9.66 2.5GeV 1.13 1.44 4.82 

- 

0.0040 - 0.002 

1.09 - 0.02 0.02 0.03 


In order to check our method, we compare our results 
with those of the ATLAS collaboration [72] and Ref. [4]. 
The comparisons are shown in Fig. gb) and Fig. gc), 
respectively. Our results are consistent with those results 
in [lElj. Figure |^c) shows the exclusion contours at 
68 % C.L. with C = 300 fb~^ for the case (b) of anomalous 
couplings. 

Conclusion and Discussion: Now equipped with the 
cut efficiency function A(mhh), we are ready to estimate 
the potential of excluding the FNNP band at the LHC. 
The SM backgrounds include 6677, CC77, bbjj, jj77, 
bbjj, ti{^ H*^), ttj, Z(—r bh)h{-^ 77), 77) and 

bbh{—r 77), etc. There are 4.72 background events af¬ 
ter all the cuts at the 14 TeV LHC with an integrated 
luminosity of 300 fb“^ [Tlj. We calculate the 95% C.L. 
exclusion bound with 


—2 In ^ —- — risj = 1.96, (16) 

where rig and nt, denotes the number of signal and back¬ 
ground events, respectively. 

Figure displays the contours of 95% C.L. exclu¬ 
sion bound from the Higgs boson pair production with 
£ = 100 fb“^ and 300 fb“^ at the 14 TeV LHC for the 
four choices of anomalous couplings. The exclusion con¬ 
tours of the cases (a) and (b) are different from those 
Rhh contours shown in Fig. The differences occur in 
the large positive C( region where BR(/i —)■ 77 ) is highly 
reduced. A large Cg is needed to reach the same rate of 
the Higgs boson pair production. 

In the case (c) the shapes of the exclusion bounds are 
very similar to those of the Rhh contours; see Fig. ic) 
and Fig.j^c). It can be understood from the fact that the 
C 2 h is not sensitive to the kinematic cuts or Br(h —r 77 ). 

In the case (d), the shape of exclusion bounds in 
Fig. I^d) is slightly different from that of the Rhh con¬ 
tours in Fig. [^d), especally in the small C 3 region. The 
difference can be understood from the cuts we imposed. 
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FIG. 5. The 95% C.L. exclusion contours in the plane of 
anomalous couplings obtained from the Higgs boson pair pro¬ 
duction with luminosity C = 100 fb“^ and 300 fb“^ at the 14 
TeV LHC. The figure indices correspond to the benchmark 
choices of anomalous couplings in Eq. 


The C 3 coupling contributes sizeably in the small mhh 
region. However, the hard px cut on bh and 77 pairs de¬ 
mands a large rrihh region. As a result, the exclusion con¬ 
tours depend mildly on C 3 . It is consistent with Ref. [55] . 

If no excess were observed in the Higgs boson pair pro¬ 
duction, then one can impose constraints on the anoma¬ 
lous couplings, especially on the FNNP band. The mini¬ 
mal luminosities to exclude the FNNP band at 95% C.L. 
in the four choices of anomalous couplings are 

(a) : £ > 1681 fb’^ (6) : £ > 235 fb”\ 

(c) : £ > 446 fb’^ (d) : £ > 186 fb”\ 

It is worth mentioning that the current allowed bands 
shown in Fig. [T] would shrink when the LHC Run-2 data 
comes out. It is possible to exclude the FNNP band with 
a luminosity smaller than the values shown above. 

Note: While this article was in finalization, the paper by 
B. Batell et al [75] appeared online and investigated the 
correlation between single Higgs and double Higgs boson 
production in supersymmetric model. 
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